Superconducting current-phase relation in Nb/Au/(001)YBa2Cu3Ox hetero junctions 
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We probe directly the superconducting current-phase relations (CPR) of Nb/Au/ (001) YBa2Cu30a: 
heterojunctions prepared in epitaxial c-axis oriented YBa2Cu30a: (YBCO) thin films. The mea- 
surements are performed by a single-junction interference experiment. Atomic Force Microscope 
observations and estimations made from I — V characteristics confirm the c-axis current transport. 
No deviation of the current-phase relation from a sinusoidal shape has been found within the ex- 
perimental accuracy. Possible reasons of the absence of the second harmonic in the CPR due to a 
second order tunneling are discussed. 



The symmetry of the order parameter in high tem- 
perature aiiperconductors has been discussed since their 
discovery.EJ Several types of phase sensitive experiment^ 
have been performed, such as |j;orner-j unction SQUIDcl 
and tricrystal ring experiment J3 The results point upon 
a dominant d-wave symmetry of the order parameter. In 
contrast, the observation of a Josephson current in c di- 
rection between a heavily twinned YBa2Cu30j; (YBCO) 
and a conventional s-wave superconductor indicates the 
presence of s companent in YBCO. Recently this conflict 
has been resolvedJj It has been shown that for YBCO 
there is mixed d and s pairing with a dominant d-wave 
component. In Ibis a case, the theoretical approach de- 
veloped in Ref.l3'B is valid and one can expect a finite 
amplitude of the second harmonic (the term I2 in the 
Josephson current, 1^ = Iism(ip) + l2Sm(2(p)) of the 
current-phase relationship (CPR) for c-axis junctions as 
well as for 45° grain boundary junctions in the ab plane. 
For the latter, the|aecond harmonic of the CPR has been 
recently observedElS On the other hand, observations 
of microwave induced steps in c-axis oriented YBCO-Pb 
tunnel junctions at multiples of hf /2e imply first-order 
currents and suggest a dominaat s-wave contribution to 
the order parameter of YBCO.Q 

In this paper we report precise measurements of CPR 
in Nb/Au/ (00 l)YBa2Cu30:, heterojunctions on YBCO 
epitaxial thin films. A pure Is = Ii sin((p) dependence 
has been observed and the reasons for the absence of the 
second harmonic in the CPR are discussed. 

Epitaxial YBCO thin films with thicknesses of 150 nm 
were obtained in situ by laser deposition on (100) 
LaAlOs substrates. The YBCO films were usually 
twinned in ab plane. The temperature of the supercon- 
ducting transition was determined by magnetic induction 
measurements as Tc = 88 90 K. The YBCO film was 



in situ covered by Au thus preventing the degradation of 
the YBCO surface layer during processing and forming an 
interface with low resistance, r = R^S — lO^^Jl-cm^, 
where Rn is the normal state resistance and 5* is the 
jimction area. Direct contacts Nb/YBCO usually have 
interfaces with resistance of the order of 10^^i7-cm^. 
Photolitography and low energy ion milling techniques 
were used to fabricate Nb/Au/YBCO junctions, with 
8 20 nm Au layer and 200 nm thick Nb films (Fig. 
1). Details of the junction fabrication were reported 
elsewhere. Il3 

The surface quality of the YBCO film is very important 
when the current transport in the c-direction is investi- 
gated. High-resolution atomic force microscope (AFM) 
studies reveal a smooth surface consisting of approxi- 
mately 200 nm long islands with maximum peak-to- valley 
distance of 3 4nm (Fig. 2). The contribution of c-axis 
and a6-plane resistance to the total resistance r can be 
estimated by using the following equation: 



r = rcTab/ [Tab + tan 7), 



(1) 



where 7 is the averaged slope angle of the islands on the 
surface, and Tc and rah are characteristic interface resis- 
tances for current transport along c-axis and in a6-plane, 
respectively. They are determined bjf the Fermi velocity 
mismatch between Au and YBCO.tj Even for a typical 
Fermi momentum anisotropy of about 3 for YBCO, one 
obtains Tc one ocder of magnitude higher than Tab for 
sharp interfaces.E3 For the typical values 7 < 2° we esti- 
mate that the tunneling area in a — 6 plane is less than 
3% of the total one. This fact indicates current fiowing 
along the c-axis of YBCOpTiifi absence of Zero Bias Con- 
ductance Peaks (ZBCE^Jlj'O which is associated with 
Andreev bound statesjlj gives additional evidence for a 
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c-axis transport in our junctions.^ 

We have measured 20 single junctions with areas in 
the range of 10 x lO^m^ to 30 x SO^m^. The I -V 
curve (see curve a, Fig. 3) shows the dependence typ- 
ical of a superconducting tunnel junction with critical 
current density Jc = 1 ^ 12A/cm^, normal resistance, 
r = 3 ^ 20/if]-cm2 and IcRn = 10 90^V {Ic is the 
junction critical current). The I — V curve for the small 
voltage scale, V < 0.3mV, is shown in the inset. Fig. 3. 
Fig. 4 presents the temperature dependence of the criti- 
cal current Ic{T) normalized to Ic at 4.2 K- The differen- 
tial resistance vs. voltage dependence Rd{V) (see curve b, 
Fig. 3) exhibits a gap-hke structure at F « 1.2mV, cor- 
responding to the Nb superconducting energy gap, Ajv6- 
This structure has a BCS-like temperature dependence 
and disappears at T « 9.1 K (see Fig. 4). At higher bias 
voltage, V > Ani,, there is a deficient current in the I—V 
curve, which, indicates the tunneling regime of junction 
conductivity.EEEj 

In the BTK model (Ref. |l^) the barrier strength is 
characterized by a dimensionless parameter Z = H/hvp, 
where vp is the quasiparticle Fermi velocity and H is the 
interface repulsive potential, responsible for an interfacial 
scattering. In our particular case Z value is character- 
ized by the ratio of the Fernu-velocities between Au and 
YBCO, VF {An) /VF (YBCO) IB Simulations of the J - ^ 
curve made on the basis of the BTK model with fitting 
parameters T = 4.2K and ANb = 1.2mV for Z = 0,1 
and 5 (see dashed curves. Fig. 3) result in a good agree- 
ment with the experimental ones for Z > 5 . Since there 
is a very small difference between / — 1^-curves for Z > 5 
in the BTK model, one can correspondingly estimate the 
lower limit for Z to be 5 and the average transparency 
D — 1/(1-1- Z^) tobe4xlO~^. The average-tijapsparency 
can be presented also as D ~ {2pcl)/ {3r)j^B'c3 where pc 
- a specific resistivity of the YBCO film in c-direction, 
I - mean free path along the c-axis of YBCO. Taking 
into account typical values pc ~ lO"'^ fl ■ cm, I ^ 1 nm, 
Rn ~ 1 S' = 100/im^, one can roughly estimate D to 
be 10^"^. This evaluation is in agreement with the com- 
parison of the experimental I — V curves with the BTK 
model given above. 

Actually, a mixed order parameter symmetry with a 
dominant d-wave component ancL a significant s-wave 
component had been confirmed.Q The YBCO is or- 
thorhombic and a small difference of the magnitude of 
the order parameter in the a- and 6-direction is expected. 
The twins in ab planes of heavily twinned YBCO can 
decrease this difference because of sign changes of the 
s-wave order parameter across twin boundaries. But be- 
cause of the disorder in the YBCO surface layer and re- 
ducing of the d-wave component across the twin bound- 



ary the d-wave component is much more suppressed than 
the s-wave one. The suppression can give rise to a thick- 
ness of a coherent_s-|jstave surface layer and non-random 
tunneling current. El'c3 

Besides, an anisotropy of the pair potential in s-wave 
superconductor (for Nb see Ref. ^) enhances the first 
harmonic component of CPR. In Ref. |2^ Rae proposed 
that this effect induced the first harmonic component of 
CPR in BizSrzCaCusOg+i (BSCCO)-Pb Josephson jimc- 
tions without assuming s-wave component in BSCCO. 
Thus, one cannot expect the Ii component to be exactly 
zero. 

CPR measurements were performed using a single- 
junction interferometer configuration in which a junction 
of interest is inserted intp-a superconducting loop with 
an inductance L sa 80pHH3 By coupling the interferom- 
eter to a tank circuit through the mutual inductance, we 
measure the impedance of the tank circuit-interferometer 
system as a function of the external magnetic flux of the 
interferometer <I>e . As shown in Ref. |2^ the CPR can be 
extracted from the following equations: 



ip^ipc- /3/((/3), 



tan a 



1 + Pf'iv) ' 



(2) 



(3) 



where (^e — 27r$e/<I>o, a is the phase angle between the 
drive current Irf and the tank voltage U at the resonant 
frequency, j3 = 27rL/c/$o is a normalized critical current 
and Q is the quality factor of the tank circuit. The details 
of this technique were reported in Ref. 27. Note, that 
from Eqsjl and ^ follow: 



tana k, k Q(3 



dm 

d(pe 



(4) 



* Theory predicts the appearance of ZBCP in the cas&.|0£ 
facets on the surface even for (100)- oriented YBCO films.EjO 



Therefore, the method outlined above is a differential 
method in respect to the CPR. This provides a high 
sensitivity of the CPR measurement. The second har- 
monic component of the CER can be detected with an 
amplitude down to lOOnA.E^ Thus, for investigated c- 
axis Nb/Au/YBCO junctions we estimate that a second 
harmonic component much lower than the first one can 
be detected, I2/I1 > 10~^. A typical CPR dependence 
measured at T = 5K is presented in Fig. 5. No devia- 
tion of the CPR from a sinusoidal 27r periodic function 
was found in any of the four samples that were investi- 
gated within the temperature range 1.6 K < T < 8.5 K. 
It means that the Josephson current is due to first ocder 
tunneling which agrees with microwave experiments.B It 
should be pointed out, that our method is more sensitive 
than in Ref. ^. 

Quite generally, 12/^1 ~ D, where D is the trans- 
parency of the tunnel barrier. If the transparency is low 
(tunneling limit) I2 is negligible compare to Ji. To be 
able to measure the second harmonic component I2, a 
tunnel junction with high transparency is required. 
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Let us discuss the difference in CPR for c-axis 
Nb/Au/YBCO and asymmetric 45° bicrystal junctions, 
investigated in Ref. |[ 

An asymmetric 45° grain boundary Josephson junction 
prepared on a bicrystal substrate satisfies such a demand 
of high transparency, as documented in Ref. ^. In that 
case, a Fourier analysis showed that /„ is non-negligible 
also for n > 3 which indicates a transparency D close to 
1 and the CPR measureipents exhibit a robust tt peri- 
odic component of Js(0).tl The formation of the Andreev 
bound states mentioned rabove even enhances the second 
harmonic component 

There are no Andreev bound states for c-axis tunnel 
junction. In the c-direction of YBCO a short coherence 
length and a Fermi velocity mismatch between low-Tc 
and high-Tc superconductors result in the transparency 
much lower compare to ah plane junction. 

From the normal-state interface resistivity (r) one can 
estimate the ratio between transparencies of the grain 
boundary [Dqb) and the c-axis (-De) tunnel junction, 
Dgb/Dc ~ 10^. Such a significant difference between 
transparencies can explain the absence of a second har- 
monic component in c-axis tunnel junction as the mag- 
nitude I2 ^ DIi is below the resolution of the method. 

The clear observation of a Josephson current with a 
pure sm{(p) CPR in a Nb/Au/YBCO tunnel junction 
can be interpreted as an indication of a dominant ^wave 
component in the order parameter of bulk YBCO.EJ The 
reason for such a conclusion, that often is given in the 
literature, is that for a pure d-wave order parameter, the 
symmetry of the c-axis tunneling between s-wave and 
d-wave superconductors leads to Ji = O.El On the other 
hand, the I2 term results from higher-order tunneling 
processes and one can neglect its weak angular depen- 
dence, i.e. I2 ^ 0. It is present also in weak links based 
on conventional s-wave superconductors, but in the tun- 
nel limit I/2/AI ^ 1- However, the CPR is 2tt periodic 
and we can exclude that current is due to higher-order 
tunneling. 

To conclude, we have found the first order tunneling 
of the Josephson current of Nb/Au/YBCO tunnel junc- 
tion due to a non-negligible s-wave component in YBCO. 
We have shown experimentally, that the second harmonic 
component I2 corresponding to the second order tunnel- 
ing, if present, is smaller than 1% of Ii. The absence of 
the second harmonic component of CPR can be caused 
by a low transparency of the c-axis Nb/Au/YBCO tunnel 
junction resulting from Fermi velocity mismatch between 
Nb and YBCO in the c-direction of YBCO and by the 
suppression of the c?-wave componcptjofpthe order param- 
eter in the surface layer of YBCO.EilH'N 
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